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The universality of fc-dependent electron-photon and electron-phonon matrix elements is 
discussed for graphene nanoribbons and carbon nanotubes. An electron undergoes a change in 
wavevector in the direction of broken translational symmetry, depending on the light polar- 
ization direction. We suggest that this phenomenon originates from a microscopic feature of 
chirality. 
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Recently, graphene nanoribbons (GNRs) have at- 
tracted attention for a variety of reasons. 1-3 ) GNRs are 
considered as unrolled carbon nanotubes (CNTs), and 
a characteristic of GNRs (CNTs) is induced by the ex- 
istence (absence) of edges. 4 ~ 6 ) GNRs as well as CNTs 
are categorized by chirality, 7 ) and armchair GNRs (A- 
GNRs) and zigzag GNRs (Z-GNRs) are known to have 
a high symmetry [see Fig. 1(a) and (b)]. Several experi- 
mental groups have attempted to clarify the optical prop- 
erties of GNRs. 8 ' 9 -' Knowing the optical selection rule 
for GNRs will be an important step in understanding 
GNRs. Theories suggest that A- and Z-GNRs exhibit dif- 
ferent polarization dependence as regards optical absorp- 
tion. 10 ~ 12 ) Namely, the selection rule for Z-GNRs pos- 
sesses a 90° rotation of polarization with respect to the 
selection rule for A-GNRs. This polarization dependence 
originates from the chirality dependent electronic wave- 
function for GNRs. 4-6 ^ However, the derivations of the 
chirality dependence described so far are complicated, 
inaccessible, and do not make it easy to grasp the essen- 
tial point of the problem. The purpose of this paper is 
to explain the selection rules for A-GNRs, Z-GNRs, and 
CNTs in a unified manner. The procedure adopted in this 
paper is not only simple but also applicable to obtaining 
selection rules even for electron-phonon interactions. 

First, we explain that the optical selection rule for A- 
GNRs is the same as that for CNTs. The electron-photon 
interaction is given by A(r, t) ■ J, where A(r, t) and J are 
a vector potential and a current operator, respectively. 
In the dipole approximation, we assume a current oper- 
ator of the form J = — (e/m)p, where — e is the electron 
charge, m is the electron mass, and p = — i/iV is the mo- 
mentum operator. This current operator does not take 
the degree of sublattice (chirality) into account. Later 
we provide a proper matrix element for the current oper- 
ator that takes account of the chirality of graphene. The 
importance of chirality will be clearly understood with 
this line of argument. 

Let us review CNTs. As a result of the cylindrical ge- 
ometry of CNTs, the wavevector (k) around the axis of 
a tube is a good quantum number, and the wavefunction 
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Fig. 1. (a) Chirality of a GNR. A large aspect ratio L^/Lz 3> 1 
(Lz/L\ 2> 1) defines an armchair (zigzag) GNR, where Lj± and 
Lz denote the lengths of the armchair and zigzag edge, respec- 
tively, (b) The chirality of a CNT is defined by the angle 8 with 
respect to the direction of the CNT axis. 



of an electron can be taken as a plane wave if we omit 
the sublattice of graphene, 
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where x and L are the circumferential coordinate and 
length, respectively. The wavevector is quantized by the 
periodic boundary condition Vfe( x + L) = ipk{x) as 
k n L = 2nn (n are integers), so that the electronic modes 
are labeled by the band index n as ip n {x). The optical 
selection rule for the nth and mth states is related to 
the possible change in wavenumber, An = m — n. The 
cylindrical geometry of CNTs gives rise to a positional 
^-dependence for the vector potential A x , 13 ^ that is, for 
a light whose polarization direction is perpendicular to 
the CNT axis, 



A x (x) = A sin ( —x 



(2) 



Then we see that An must be ±1 in order to have a non- 
zero matrix element of — (e/m) § A x (x)Tp^ n (x)p x ip n (x)dx . 
This condition for wavenumber An = ±1 is the selection 
rule of CNTs with a perpendicularly polarized light. 13-15 ) 
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With GNRs, a plane wave is reflected at the edge. 
Consequently a standing wave is formed by the super- 
position of two plane waves propagating in opposite di- 
rections. Suppose that an incident wave has the wavevec- 
tor k in the direction perpendicular to the edge, then the 
wavevector of the reflected wave is given by —A: as a result 
of momentum conservation. There are two possible types 
of superpositions: ipk{x)±ip-k(x). For A-GNRs, the elec- 
tron wavefunction should vanish at the edge (x = 0), and 
only the antisymmetric combination, 



(3) 



<Pk{x) = y £ sin(fca:), 



is selected. The boundary condition should be also im- 
posed at x — L as (flk{L) = 0, which gives k n — irn/L 
(n are positive integers). The standing waves are la- 
beled by the band index n as tp n (x). The flat geom- 
etry of GNRs results in a constant vector potential 
A x (x) = A, and the matrix elements are given by 
— (e/m)A J Q ip m (x)p x (p n (x)dx. This integral results in 



sin(fc m a;) coa(k n x)dx 
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m — n m + n 



n G even 



m — n € odd. 



(4) 



Because the Fermi wavevector satisfies fcp = 47r/3a (a 
denotes a lattice constant), the values of m and n are se- 
lected so that they are large enough for (m + n) -1 to be 
negligible compared with (m — n) -1 . Thus, we approxi- 
mate ji J Q sm(k m x) cos(k n x)dx — (2/-K)An~ 1 . Further- 
more, the transition amplitudes for An = ±3, ±5, . . . are 
suppressed by the factor of An -1 . The selection rule of 
A-GNRs is mainly given by An = ±1, which is coinci- 
dent with the selection rule of CNTs as discussed above. 

It is meaningful to consider the reason for the sim- 
ilarity. For CNTs, the change in electron wavenumber 
(An = ±1) is brought about by the macroscopic topol- 
ogy of the cylinder. A light behaves as a standing wave 
on a cylindrical surface [eq. (2)], while the electrons are 
plane waves [eq. (1)]. For GNRs, the electrons are stand- 
ing waves [eq. (3)], while the light is a uniform plane 
wave (i.e., zero- wavenumber mode). An essential factor 
in obtaining An = ±1 for A-GNRs is that the momen- 
tum operator p x alters sm(kx) into cos(fcx). Since the 
momentum operator originates from a lattice spacing as 
ip n (x + a/2) — ip n (x — a/2) w a(d / dx)<p n (x) , we see that 
a microscopic lattice topology is essential for a change 
in the electron wavenumber. Contrastingly, a plane wave 
satisfies a(d / dx)il> n (x) oc tl> n {x), and the orthogonality 
condition tells us that the corresponding current cannot 
change the electron wavenumber when A x (x) is constant. 
Therefore, as long as standing waves are formed in GNRs, 
we expect the wavnumber to change in an optical tran- 
sition. A change in the wavenumber will be a universal 
phenomenon for GNRs, as we see in the following. 

To investigate the polarization dependence of the se- 
lection rule, we now take account of the local arrange- 
ment of the carbon-carbon bonds on an atomic scale [see 




Fig. 2. (a) The bond vectors, ri = a cc e y , r-2 = a cc (— -^-e x — 

^e y ), and r3 = a cc (-^-e x — ^e^), where a cc denotes the bond 
length, (b) The band index s in the Dirac cone, and the relative 
amplitude between two sublattices is determined by s and ©, 
respectively. These variables correspond to the amplitudes at 
two sublattices. 



Fig. 2(a)]. Since graphene's hexagonal unit cell consists 
of two atoms (A and B), the wavefunction has two com- 
ponents 



(r A ) = 



*tfc BS ( r A) 

*Lfc,a(rA + ri; 



(5) 



where s is the band index (+1 for the conduction band, 
— 1 for the valence band), ta = (x,y) denotes the posi- 
tion of an A-atom, r a (a = 1,2,3) represents the bond 
vectors [see Fig. 2(a)], and 6 is the polar angle in (k x ,k y ) 
plane defined with respect to the Dirac point at (fcF,0), 
as shown in Fig. 2(b). The current also consists of two 
components: J A (rA) and J B (re), where J A (rA) [J B (re)] 
represents the electron's flow into an A-atom [B-atom] 
at ta [tb] from the nearest-neighbor B-atoms [A-atoms] . 
These components contribute to a local dipole moment 
and are essential to the optical transition. 16 ) An opti- 
cal transition in each sub-atom does not induce a local 
dipole moment and can be omitted from the analysis. 

Suppose that an incident light transfers an electron 
from the B-sites denoted by open circles in Fig. 2(a), 
to the central A-site (solid circle). The corresponding 
current component is given by 

J Lfe; s ',fc„fe yS ( r A) - ie^ mk , yS ,(r A )^x 

I E r a ^ feyS (r A + r a )l, (6) 

la=l,2,3 J 

where 7 denotes the nearest-neighbor hopping inte- 
gral. By putting eq. (5) into eq. (6), we obtain 

/ J Lk' y s',k n k y s(*) d y = «fctfc»fe <e( * m,fcy) jmn(a:). where 
itn(x) =ie^ m {x) x \r lVn {x) +r 2 e- lfe «( 3a -/ 2 V„(x- |) 

+ T3e -ik y (3a cc /2) ipn{x+ a ) y (?) 

Because e -*M 3a cc/2) j g approximately equal to 1 for the 
states near the Dirac point, we set e -*M 3a cc/2) = 1 m 
eq. (7). 12 ) Thus, k y does not appear explicitly in the 
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following analysis for A-GNR, but it is implicitly taken 
into account through Q(k n ,k y ). Here let us define the 
amplitude for an A-atom of GNR [aj] mn = / Ai(x)et ■ 
3mn( x )dx (i = x,y), which are given by 

r-L 



(8) 



\oLx\mn = iAevFy / sm(k m x) cos(k n x)dx, 
L Jo 

2 f L 

[a y ]mn — iAev-p— J sm(k m x)sm(k n x)dx, 

where v F = 3ja cc /2H. It is straightforward to show 
[a x ]mn = iev F A(2/n)An^ 1 by using eq. (4), and 
[«y]mn = iev-pA5 mn owing to the orthogonality condi- 
tion. Similarly, the local electronic current component at 
a B-atom is written as 

I r Q < feHs(rB _ ra) l (9) 

lo=l,2,3 J 

where re denotes the position of a B-atom. We de- 
fine the current and amplitude for a B-atom as 
fJLk>s>, knkyS (rB)dy = S k , ky ^e-^^jta(x) and 
[ft]mn = / Ai(x)&i ■ j^ ln (x)dx, respectively. We obtain 

2 f L 

[Px]mn = iAev F y / sm(k m x) cos(k n x)dx, 

L ™ 

2 f 

[f3y} mn = -iAev F — J sm(k m x) s'm(k n x)dx. 

The optical transition amplitude for A-GNR is con- 
structed from the sum of the current components, 
II J Lk yS \k n k v s( r A)dxdy and JJ Jf mkyS , knkyS (r B )dxdy, 



as 



Mi(k m kyS , k n kyS 
1 



By comparing eq. (8) with eq. (10), we see that [f3] and 
[a] are related to each other via 



IP: 



x\mn — l/^xjmm 



tt 



y\mn- 



(12) 



Then, the matrix element of eq. (11) may be written in 
a more compact form (in evpA units) as 

M x (k m k yS ',k n k vS ) = -J- (se m ^^ + s'e- i@ ^ 
My(k m kyS ,k n kys) = 5 mn — ^se 



i&(k m ,k v ) _ s > e -i&(k n ,k v ) 



(13) 

Equation (13) is the optical matrix elements for A-GNRs. 
This result was used to explain the optical absorption 
spectra in ref. 12. 

The matrix elements for Z-GNRs are obtained by re- 
peating a similar calculation to that given above. A 
standing wave for Z-GNRs is formed by the superpo- 
sition of two waves propagating in opposite y-dircctions. 
Suppose that an incident wave has a wavevector k y in a 
direction perpendicular to the edge, then the wavevector 



of the reflected wave is given by — k y as a result of mo- 
mentum conservation. The possible superpositions are 

1 / e -i&(k x ,k v )\ l / e -i0(k x -k v )^ 

(14) 



Here, we assume that the zigzag edges appear at y = 
and y — i, and consider a case where the zigzag edge at 
y = (y = L) consists of A-atoms (B-atoms) . Then, the 
electron wavefunction at the B-atoms should vanish at 
the edge (y = 0), and the antisymmetric combination is 
selected. Thus, the standing wave for Z-GNRs is written 
as 

k s (x,y) = ^ k (x)4= 

k x k v s\ ,yj vk x \ ssm(k y y) 

■ (15) 

The boundary condition should also be imposed at y = L 
as ^ k k s (x,L) — 0, by which k y is quantized accord- 
ing to k n L — Q(k x ,k n ) — nir. Using the wavefunction 
^k x k„s(x, y), we obtain the optical matrix element, 12 ) 



^x{k x k m s , k x k n s) 



=- (scose, 



s' cos 0„) [a y 



(16) 



- - (s sin 9 m [a x ] nm + s' sin 6„ [a x 
My{k x k m s , k x k n s} 

= i (s cos 6 m - s' cos 6„) [atylmn 

- i (ssmQ m [a x ] nm - s'sme n [a x 



where we abbreviate @(k x ,k n ) as 9„ and use the fact 
that k x w kp to obtain the right-hand side. For an inter- 
band transition (s = — 1 and s' = 1), we see that the 
wavenumber selection rule is given by An = ±1, ±3, . . . 
(An = 0) when the polarization of an incident light is 
set parallel (perpendicular) to the zigzag edge, which is a 
90° rotation of polarization with respect to the selection 
rule for A-GNRs. For an intra-band transition (s = s'), 
the selection rule for wavenumber exhibits a 90° rotation 
of polarization with respect to the An-selection rule for 
an inter-band transition. This additional change in the 
polarization dependence is a characteristic feature spe- 
cific to Z-GNRs, not shown in the An-selection rule for 
A-GNRs. 

In Table I, we show our results for the electron-photon 
matrix elements for GNRs, where we pay special at- 
tention to the change in wavenumber An. By replacing 
ip n (x) in eq. (5) with tp n (x), we see that the polarization 
dependence of the An-selection rule for CNTs is irrele- 
vant to the chirality. The existence of the edge and the 
resultant formation of the chirality-dependent standing 
wave are important to obtaining the chirality-dependent 
An-selection rule. 

Finally, we briefly mention that this method is ap- 
plicable to an electron-phonon interaction. A relative 
displacement vector for an optical phonon mode u(r, t) 
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Table I. Polarization dependence of An for the possible inter- 
band optical transition, "para" ("perp") means that the polar- 
ization is parallel (perpendicular) to the edge (axis) of a GNR 
(CNT). For GNRs, An induced by an electron-phonon interac- 
tion follows the same rule, where "para" ("perp") means the 
direction of the vibration. 





A-GNRs 


Z-GNRs 


CNTs 


para 


An = 


An = ±1,±3, . . . 


An = 


perp 


An = ±1,±3,.. . 


An = 


An = ±1 



undergo a change in wavenumber through an optical, 
inter-band transition when the polarization of an inci- 
dent light is perpendicular (parallel) to the armchair 
(zigzag) edge. The origin of the shift in wavenumber 
is attributed to the microscopic topology of the lattice, 
in that the standing wave has different probability am- 
plitudes at adjacent sites in a chirality dependent man- 
ner. Thus, a change in wavenumber is a universal phe- 
nomenon in GNRs. The universality is also recognized 
for electron-phonon interactions. 

K.S acknowledges a MEXT Grant (No. 23310083). R. 
S. acknowledges a MEXT Grant (No. 20241023). 



(= ua - ub) changes the hopping integral from 7 to 
7 + gu(r,t) • (r a /a cc ), where g is a coupling constant, 
and induces a deformation current. The current induced 
at an A-atom by the lattice deformation u(r) is given 
as J^' s ' ks( r A-) by replacing —67 with g/a cc . The current 
induced at a B-atom is given by — Jk' S ',ks( r B)> where 
the extra minus sign originates from the fact that u(r, t) 
changes its sign when we replace A and B-atoms: the in- 
teraction must be symmetric with respect to the change 
of A and B as (J A - J B ) • (ua - ub). Therefore, the 
transition amplitude is constructed from the difference 
between the current components, 

J f |Jk< s <, ks (rA) - JkV,k s (r B )} dxdy. (17) 

Since the minus sign is taken into account as s' — > — s' in 
Eqs. (13) and (16), we conclude that the electron-phonon 
matrix element for an intra-band (inter-band) transition 
is the same as the electron-photon matrix element for an 
inter-band (intra-band) transition, except for a numeri- 
cal factor. Thus, the selection rule for the wavenumber 
of a phonon mode whose vibrational direction is either 
parallel or perpendicular to the edge allows a non-zero 
shift An = odd in the wavenumber of an electron (see 
Table I). 

In conclusion, the electronic standing waves in GNRs 
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